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Global quantum discord (GQD), proposed by Rulli and Sarandy [Phys. Rev. A 84, 042109 
(2011)], is a generalization of quantum discord to multipartite states. In this paper, we provide 
an equivalent expression for GQD, and obtain the analytical expressions of GQD for two classes of 
multi-qubit states. The phenomena of sudden transition and freeze of GQD are also discussed. 
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I. INTRODUCTION 

Quantifying the multipartite quantum correlations is a 
very challenging and still largely open question [1-5] . For 
bipartite case, entanglement and quantum discord have 
been widely accepted as two fundamental tools to quan- 
tify quantum correlations [6, 7], and quantum discord 
captures more quantum correlations than entanglement 
in the sense that a separable state may have nonzero 
quantum discord. Generalizations of bipartite quantum 
discord to multipartite states have been considered in dif- 
ferent ways [8] . In [5] , Rulli and Sarandy proposed a mea- 
sure for multipartite quantum correlations, called global 
quantum discord (GQD), which can be seen as a gener- 
alization of bipartite quantum discord [9, 10] to multi- 
partite states. GQD is always nonnegative and its use 
is illustrated by the Werner-GHZ state and the Ashkin- 
Teller model [5]. 

In this paper, we provide an equivalent expression for 
GQD, and give an interpretation of GQD (Sec. III). We 
derive the analytical expressions of GQD for two classes 
of multi-qubit states (Sec. IV), these results generalize the 
earlier results [9, 11]. The phenomena of sudden transi- 
tion and freeze of GQD are also discussed (Sec.V). For 
clarity of reading, we first recall the definition of GQD 
proposed in [5] (Sec. II). 



II. GLOBAL QUANTUM DISCORD (GQD) 



We briefly review the definition of GQD proposed in 
[5]. 

Consider two systems Ay and A 2 (each of them is of 
finite dimension), the symmetric quantum discord of a 
state pa x a 2 °f the- composite systems A\A 2 is 

D (PA 1 A 2 ) = mm[I(p AlA2 ) - I(^ A1 A 2 (PA 1 A 2 ))}- (1) 



In Eq.(l), 

I (PA 1 A 2 ) = S(p Al ) + S(pA 2 ) ~ S{p Al A 2 ), 



(2) 



is the mutual information of p Al A 2 , min is taken over all 
locally projective measurements performing on AB, <£>(.) 
denotes a locally projective measurement performing on 
the system (•), S(-) is the Von Neumann entropy, and 
PAi, PA 2 are reduced states of p Ai A 2 - 

D(pa 1 a 2 ) is a natural extension of the original defini- 
tion of quantum discord which defined over all projective 
measurements performing only on A\ or A2 [9, 10]. 

Since the mutual information I(pa 1 a 2 ) can be ex- 
pressed by the relative entropy 

I (PA 1 A 2 ) = S{p Al A 2 \\PA 1 ®PA 2 ), (3) 

hence, Eq.(l) can also be recasted as 

D (PA t A 2 ) =min{S{p Al A 2 \\pA 1 ® PA 2 ) 
S(® Al a 2 (pa 1 a 2 )\\$ > a 1 (pa, ) ® $a 2 [pA 2 ))] ■ (4) 

Note that the relative entropy of state p with respect to 
state a (p and a lie on the same Hilbert space) is defined 
as 

S (fi\\(r) = tr(p\og 2 p) - tr(p\og 2 a). (5) 
Further, Eq.(l) can also be rewritten as 

D (PAxA 2 ) = wn[S(p Al A 2 ll'^A 1 A 2 (pA 1 A 2 )) 
-S( PAl ll$ Al ( PAl )) - S(PA 2 I^A 2 (PA 1 ))]- (6) 

The definition of GQD is a generalization of bipartite 
symmetric quantum discord. Consider N (2 < N < 00) 
systems Ai, A% , ... , An (each of them is of finite di- 
mension), the GQD of state pa 1 a 2 ...a n on the composite 
system AiA 2 ...An is defined as [5] 



D (PA 1 A 2 

= min[S , (p j4l A 2 

N 



r) 



A X A 2 



,(PA lJ 



,)) 



-J2S(pA J \\^A J (PA j ))}- 



(7) 
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It has been proved that D(pa 1 a 2 ...a n ) > for any state 
Pa x a 2 ...a n [5]. Also, it is easy to see that D(p Al A 2 ...A N ) 
keeps invariant under any locally unitary transformation. 
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III. AN EQUIVALENT EXPRESSION FOR 
GLOBAL QUANTUM DISCORD 

In this section, we provide an equivalent expression 
for GQD. We first state two mathematical facts as the 
lemmas below. 

Lemma 1. For any square matrix (with finite dimen- 
sion) A, let A be the matrix whose diagonal elements are 
the same with A, and other elements are zero. B and B 
are defined similarly. Then 

tr(AB) = tr(AB), (8) 
tr(Af(B)) = tr(Af(B)), (9) 

where /() is any function. 

Lemma 2. Let PAiA 2 ...a n be a state on Hilbert space 
H12...N, PAn PA 2 , ■■■ > PA N be the reduced states of 
Pa 1 a 2 ...a n on Hilbert spaces Hi, H2, ... , Hpj, respec- 
tively. Suppose a An ffA 2 , va n are states on Hi, Hi, 
... , Hn, respectively. Then it holds that 

tr[p Al A 2 ...A N \og 2 {a Al ® <JA 2 <2> ••• <2> <JA N )} 

N 

= ^2tr i [p Ai log 2 a Ai ]. (10) 

i=l 

Proof. We only prove the case of N — 2, the proof of 
N > 2 is similar. When N = 2, we need to prove 

tr[p AlA2 \og 2 (a Al ®<7A 2 )\ 
= tn [p Al log 2 cr Al }+ tr 2 [pa 2 log 2 <7a 2 ] ■ (11) 

It is known that p Al A 2 can be written as [12] 

PMM = c jPij ® P2j, (12) 
j 

where {cj}j are real numbers, py, p 2j are all Hermite 
matrices. For o~ Al , o~a 2 , there exist unitary matrices Ui 
and U2 such that D\ = \J\oaJJi , D 2 = U 2 ga 2 U 2 arc 
all diagonal, where + denotes adjoint. Note that 

log 2 (L>i ® D 2 ) = (log 2 Di) ® I 2 + h ® (log 2 D 2 ), (13) 

where 7i, 7 2 arc the identity operators on Hi , H 2 , re- 
spectively. Then Eq.(ll) can be directly verified. □ 

With the help of Lemma 1 and Lemma 2, we can get 
an equivalent expression for GQD defined by Eq.(7). 

Theorem 1. The GQD of a state pa 1 a 2 ...a n defined by 
Eq.(7) can also be expressed as 

D(pa 1 a 2 ...a n ) 

= mm[I(p Al A 2 ...A N ) - I($A 1 A 2 ...A N (PA 1 A 2 ...A N ))}, (14) 

where, the mutual information 

N 

I(PA 1 A 2 ...A N ) = Y S (PAi) - S{p Al A 2 ...A N )- (15) 

t=l 



Proof. From lemma 1, we have 

S{ PAj \\$> Aj {p Aj )) = -S{p Aj ) + S(* Aj {p Aj )). (16) 
Note that 

N 

I{PA 1 A 2 ...A N ) = Y S (P A ^ _ S(p Al A 2 ...A N ) 
i=l 

= S(p Al A 2 ...A N \\PA 1 ®PA 2 O PA N )- (17) 

Together with Lemma 2, we can easily prove Theorem 1. 
□ 

From Theorem 1, we see that, GQD of a state is just 
the minimal loss of mutual information over all locally 
projective measurements. This interpretation of GQD 
is consistency with the symmetric quantum discord in 
Eq.(l), as well as the original definition of quantum dis- 
cord for bipartite states. 

For a special case, we consider a state p Al A 2 ...A N whose 
reduced states p Al , pa 2 , ■■■ , Pa n are all proportional to 
identity operator. In such case, the GQD of pa x a 2 ...a n 
can be remarkably simplified. We state it as Theorem 2, 
its proof is easy. 

Theorem 2. An TV-partite state pa 1 a 2 ...a n , if its re- 
duced states pAn Pa 2 , ■■■ , PA N are all proportional to 
identity operator, then the GQD of p Al A 2 ...A N can be 
expressed as 

D(pA 1 A 2 ...A N ) = ~S(pA 1 A 2 ...A N ) 

+ mmS(^> Al A 2 ...A N {pA 1 A 2 ...A N ))- (18) 



IV. TWO CLASSES OF iV-QUBIT STATES 
WHICH ALLOW ANALYTICAL EXPRESSIONS 
OF GQD 

We consider two classes of iV-qubit states which allow 
analytical expressions of GQD. We first recall two math- 
ematical facts. 

Lemma 3. [13] Group homomorphism of U(2) to 
SO(3). 

For any two-dimensional unitary matrix U, there ex- 
ists a unique real three-dimensional orthogonal matrix 
R with determinant 1, such that for any real three- 
dimensional vector f= {r x ,r y ,r z ) , it holds that 

Uf-aU+ = R(f)-a. (19) 

Conversely, For any real three-dimensional orthogonal 
matrix R with determinant 1, there exists (not unique) 
a two-dimensional unitary matrix U, fulfills Eq.(19). 

In Eq.(19), f-a = r x a x +r y a y + r z a z , a = {o- x ,a y ,a z } 
are Pauli matrices; R(r) = (i?f*)* is a real three- 
dimensional vector, here t denotes matrix transpose. 

Lemma 4- [14] Monotonicity of entropy function under 
majorization relation. 

For given {pi,p 2 , -,p n }, {Qi,Q2, -,q n }, satisfy 

1 > Pi > P2 > ■■■ > Pn > 0, Yn=iPi = !. 



3 



1 > qi > Q2 > > In > 0, £" =1 % = !• If 

fe fc 

<^2<lt,k = 1,2,..., n, 



i=i i=i 



then 



(20) 



(21) 



i=i 



Theorem 3 and Theorem 4 below provide two classes 
of multi-qubit states which allow analytical expressions 
of GQD. 

Theorem 3. For N-qubit (N > 2) Werner-GHZ state 



Y 1 



p = (! -m)-^ + mI^><>| 5 



(22) 



the GQD of p is 



l-,i 



^(P) = (-j^r + M) ^(-^ + m) + -gjvr 1°8 2 "^jv" 



o, 1 - A* , Ms , ,1 - M . AS 



(23) 



Where, I is 2 x 2 identity operator, \x G [0, 1], |V>) is the 
N-qubit GHZ state 



|V> = (|00...0) + |ll...l>)/>/2. 



(24) 



Proof. A projective measurement on single qubit can 
be expressed as 

n = i(/ + II.a), U^^I-U-a), (25) 

where, II = (a, /3, 7) is a real vector with unit length. 
It is easy to check that 

U(a x ) = U a x U +Il 1 a x U 1 =an-a, (26) 
n(o- v ) = Iloo- y Tlo+Il 1 o- y Il 1 =f3ii-a, (27) 

u(a z ) = n CT z n + nx^n! = 7 n • s. (28) 

Now for any locally projective measurement performing 
on A^-qubit state p, we label it by IIi = (ai, /3i, 71), 

n 2 = (02,^2,72), •■• j IIjv = (ajv,y8jv,7Jv), each of them 
is a real vector with unit length. It is easy to verify that 
for state p in Eq.(22), whose reduced states pA t , Pa 2 ? ••• > 
PAn are all proportional to identity operator. So, we can 
calculate -D(p) according to Theorem 2. We then need 
to calculate S(p) and S($ Ai a 2 ...a n (pa 1 a 2 ...a n ))- 
The eigenvalues of p can be easily found, that is 

5pec(p) = {^+ M ,^,^r,...,^r}. (29) 
can then be directly calculated by Eq.(29). 



Let < J > AiA 2 ...Aiv be the locally projective measurement 
{n»}£Li, using Eqs. (26-28) and the facts 



^(|0)(0|® Ar + |0)(l|® Ar 



+ |l)(0| 8A, + |l)(lO, (30) 



|o)(o| = ^, |i)(i|= 7 ^ 



|o)(i| = ^±^, |i><0| 



2 ' 

C7l - i<7 2 



2 ' 1 " 1 2 
we get 

$ Al A 2 ...A N (\il>)(i>\) 
1 w / + 7jIL. • <? ^ ctj + 



(31) 
(32) 



+ ® i= i o n * ' CT + ®i=l _ 



(33) 



2 * 1=1 2 

For given {I?i}-^ l5 according to Lemma 3, we can 
always perform a locally unitary transformation on 
$AiA2...A N (p) (S($a 1 a 2 ...a n (p)) keeps invariant) such 
that LTj • a is transformed to cr^ for all i, then Eq.(33) 
becomes 



2 L ~ l=1 2 



N a t + i/3i 



N ai ipi , , 
<T : »i=i ^ o z \. (34) 



1=1 2 1=1 2 

Eq.(34) is a diagonal matrix with 2 N diagonal elements 



1=1 



n o i + (-l)^ + ^ i -(-l)^ L (3g) 



where, ni, n 2 , njy € {0, 1}. Consequently, 
$AiA 2 ...a n (m) have 2^ eigenvalues 



2 JV 



1=1 



i=l 



where, m, n 2 , G {0, 1}. 

Note that in Eq.(36), the eigenvalues corresponding to 
{ni,n 2 , tin} and {1 — m, 1 — n 2 , 1 — ^at} are equal. 
Together with Lemma 4, we then assert that when 7^ = 1 
for all i, S($AiA 2 ...A N (p)) achieves its minimum, since 
when 7, = 1 for all i, $^1^2... An (p) have 2 N eigenvalues 

f l-M , M 1 - A* A* 1 ~ M 1 ~ M 1 ~ M i /o 7 n 
2 W 2' 2^ 2' 2 W ' 2 N oAr 



2 Af 



So, by Eq.(37) and Eq.(29), Theorem 3 can be proved. 
□ 
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We consider the behavior of Eq.(23) with large number 
N. Let N — » oo, we get Corollary 1 below. 

Corollary 1. When N — > oo, Eq.(23) approximates 



(38) 



We remark that, Eq.(38) is somehow a good approxi- 
mation of Eq.(23). In fact, when N = 10, the deviation is 
less than 10 -2 ; N — 14, the deviation is less than 10~ 3 ; 
N = 17, the deviation is less than 10~ 4 . Fig.l shows 
Eq.(23) for N — 2, 3, 5, oo, respectively. 

We also remark that, when N — 3, Eq.(23) is the same 
to the result obtained in [5]. When N = 2, Eq.(23) re- 
turns the result of Werner state considered in [9] . 






N=oo 




N=5 




N=3 




N=2 



when N is even, 

4 

fl (p) = -l-5^A i log a A 3 -, (45) 

3=1 

Ax = [l + C3 + ci + (-l) Ar/2 c 2 ]/4, (46) 

A 2 = [l + c 3 -c 1 -(-l) Ar / 2 c 2 ]/4, (47) 

A 3 = [l- C3 + c 1 -(-l) Ar / 2 c 2 ]/4, (48) 

A4 = [l-c 3 - Cl + (-l) JV / 2 c 2 ]/4. (49) 

In Eq.(39), / is the 2x2 identity operator, {c\, c 2 , c 3 } are 
real numbers constrained by d G [0,1] (when TV is odd) 
or Ai, A 2 , A 3 , A4 e [0, 1] (when N is even). 

Proof. Obviously, for state p in Eq.(39), its reduced 
states pa 1 , pa 2 , ■■■ , PA N are all proportional to iden- 
tity operator. Thus, we can calculate D(p) according 
to Theorem 2. We then need to calculate S(p) and 
S($ Ai a 2 ...a n {pa 1 a 2 ...a n ))- It is easy to see that state p 
has nonzero elements only on the principle diagonal and 
the antidiagonal, so is the matrix (p — xI® N ). By the 
Laplace theorem in linear algebra, M — det(p — xI® N ) 
can be expanded as the multiplications of 2 Ar_1 determi- 
nants of 2 x 2 matrices, 



0.0 0.2 0.4 0.6 0.8 1.0 



M = det(p - xl 

2 N-1 

= Jl det 



M 



j j 



M. 



l2 N -j,j M2 N -j,2 N -j 



(50) 
(51) 



FIG. 1: GQD of Werner-GHZ state in Eq.(22) when N = 2, 
3, 5, 00, respectively. 



We next consider another class of A^-qubit states. 
Theorem 4- For A^-qubit state 



where 



Al 



M 



M . . — 1 ( 1 \n 1 +n 2 + ...+n N _£3_ _ 



£1 

2 N 

Cl 



N+ ni +n 2 + ...+n N _£3_ 

2 N 
c 2 

2 N '' 



(52) 
x, (53) 



i ;N / 1 \N+m+n 2 + ...+n N C 2 /rr\ 

i,2«-_j - ^ + 1 (-1J (55) 



the GQD of p is 



In Eq.(40), 



D(p) = f(p)-g(p). 



g(p) — i»^_> 



<* = V c ? + c i + "' 



'3' 



(40) 



f{p) = — log 2 — g-loga-g-, (41) 

c = mai{|ci|, |c 2 |, |c 3 |}; (42) 

when A" is odd, 



1 + d. 1 + d l-d, l-d , JO , 
^- l0g2 ^ 2- l0g2 ^-' (43) 



(44) 



{n 1 ,n2,—,n N } € {0,1}. 

Then by direct calculation of det(p — xI® N ) = 0, we 
get the eigenvalues of p. That is, when A" is even, the 
eigenvalues of p are 



r 1 + C3 ± Cl + (-1)^2 l-c 3 ± Ci 



(-l) N/2 c 2 



2 N 



2 N 



' 2 W 



2 N 



h 



(56) 



each of them possesses multiplicity 2 N 2 . When A~ is 
odd, the eigenvalues of p are 



{^(li^ + cl+c 2 )}, 



(57) 



each of them possesses multiplicity 2 W 1 . 

S'(p) can then be directly calculated by Eqs. (56-57). 
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We now calculate &a 1 a 2 ...a n {p)- After the measure- 
ment {Ui — (oti, /8»j7»)}£Li as labelled in the proof of 
Theorem 3, from Eqs. (26-28), we have 

1 N 

*a iA2 ...a n (p) = ^[i® N + (aHm 

1 i=l 

N N 

+c 2 l[p i + c 3 l['y i )^ =1 (U i -a)}. (58) 

i=l i=l 

From Lemma 3, for any {IL}^ we can always per- 
form a locally unitary transformation on &a 1 a 2 ...a n (p) 
(S($a 1 a 2 ...a n (p)) keeps invariant) such that flj • u is 
transformed to a z for all i, then Eq.(58) becomes 

N N N 

w ii® n + (ci n ^ + <* n & + <* n ^ w ] . (59) 

«— 1 i— 1 i— 1 

Eq.(59) is a diagonal matrix, so its eigenvalues are 

N N N 

{^r [1 ± (ci n ^ + c2 n ^ + n ( 6 °) 

z— 1 i—1 i—1 

each of them has multiplicity 2 Ar_1 . 

From Lemma 4, we find that minimizing 
S(&a 1 a 2 ...a n {p)) is equivalent to maximizing 

N N N 

|ciJJa i + c 2 JJ/3i + C3jj7i| ) (61) 

i—1 i—1 i—1 

over all possible {n»}£Li- 

Suppose {ftj}^ 1 are given, then Eq.(61) can be writ- 
ten as 

N-l N-l N-l 

n^n^ ( 62 ) 

i—1 i—1 i—1 

So over all possible Un = (ajv, Av, 7jv) with + + 
7^ = 1, Eq.(62) achieves the maximum 

^n^+^n^+^n^ 2 ) 1 / 2 . (63) 

i—1 i—1 i—1 

Suppose {fli}^ 2 are given, because af v _ 1 + + 
■y%_ 1 = 1, so the maximum of Eq.(63) over all possible 
IIjv-i = (ajv-i,^jv-i,7Jv) is 

(ma X {c\ jf al c\ jf ft, c\ jf^}) 1 /*. (64) 

i—1 i—1 i—1 



The maximum of Eq.(64) over all possible {ni}^ 2 ap- 
parently is 

c = max{\d\, |c 2 |, |c 3 |}. (65) 

By Eqs. (56, 57, 60, 65), we then complete this proof. □ 

We remark that, when N = 2, Theorem 4 recovers 
the result of 2-qubit Bell-diagonal state which was first 
obtained in [11]. 

We know the original quantum discord may manifest 
the phenomena of sudden transition and freeze [15-17]. 
With the analytical result of Theorem 4, we assert that 
GQD can also manifest such interesting phenomena. We 
make this assertion more clear by giving an example. Let 
N qubits be in the state Eq.(39), and let any qubit un- 
dergo a local phase damping 

E = y/1- p/2 I, E x = ^j2o z . (66) 
After this channel, the state p in Eq.(39) becomes 

Pip) = ±{I* N + Clip)*®" + c 2 (p)af N + c 3 {p)af N ), 

(67) 

where 

ci(p) = ci(l-p), c 2 (p) = c 2 (l -p), c 3 (p)=c 3 . (68) 

So, D{p(p)) can be calculated by Theorem 4. Therefore, 
similar to the bipartite case discussed in [15, 16], it can 
be found that, the sudden transition occurs if and only if 

0< |c 3 | <max{|ci|,|c 2 |}, (69) 

and freezing GQD may occur when N is even. 

V. CONCLUSION 

In summary, we provided an equivalent expression of 
global quantum discord (GQD). From this equivalent ex- 
pression, we gave an interpretation of GQD as the mini- 
mal loss of mutual information over all locally projective 
measurements. This interpretation is consistency with 
the original quantum discord. We obtained the analyt- 
ical expressions of GQD for two classes of multi-qubit 
states, each of them possesses high symmetry. By the 
analytical expressions of these states, we discussed some 
behaviors of GQD, including the asymptotic behavior of 
GQD when N tends infinity, and the phenomena of sud- 
den transition and freeze of GQD. 

This work was supported by the Fundamental Re- 
search Funds for the Central Universities of China (Grant 
No.2010scu23002). The author thanks Qing Hou and Bo 
You for helpful discussions. 
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